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ABSTRACT 

Current observations suggest that the universe was reionized sometime before 
z ~ 6. One way to observe this epoch of the universe is through the Near Infrared 
Background (NIRB), which contains information about galaxies which may be 
too faint to be observed individually. We calculate the angular power spectrum 
(Ci) of the NIRB fluctuations caused by the distribution of these galaxies. As- 
suming a complete subtraction of any post-reionization component, Ci will be 
dominated by galaxies responsible for completing reionization (e.g., z ~ 6). The 
shape of C\ at high / is sensitive to the amount of non-linear bias of dark mat- 
ter halos hosting galaxies. As the non-linear bias depends on the mass of these 
halos, we can use the shape of C% to infer typical masses of dark matter halos 
responsible for completing reionization. We extend our previous study by using 
a higher- resolution iV-body simulation, which can resolve halos down to 10 s M . 
We also include improved radiative transfer, which allows for the suppression of 
star formation in small-mass halos due to photo-ionization heating. As the non- 
linear bias enhances the dark-matter-halo power spectrum on small scales, we 
find that C\ is steeper for the case with a complete suppression of small sources 
or partial suppression of star formation in small halos (the minimum galaxy mass 
is M min = 10 9 M Q in ionized regions and M min = 10 8 M in neutral regions) than 
the case in which these small halos were unsuppressed. In all cases, we do not 
see a turn-over toward high I in the shape of l 2 Ci. 



Subject headings: cosmology: theory - - diffuse radiation - - galaxies: high- 
redshift — infrared: galaxies 
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INTRODUCTION 



Probing the beginnings of star and galaxy formation in the early universe is one of the 
goals of modern cosmology. These high redshift stars provided a wealth of ionizing photons, 
which caused all of the hydrogen in the universe to be ionized, a process called reionization. 
Therefore, understanding high-redshift star formation is closely coupled with our understand- 
ing of reionization history. Recent measurements of the polarization of the Cosmic Microwave 
Background have shown that reionization started early an d was extended in time, with an 



equivalent i nstantaneous rei o nization redshift of z ~ 11 (IKogut et al 



2003 . 120071 : IPage et al.l 120071 iDunklev et al.ll2009l : Komatsu et all 120091 ). Therefore, we ex 



20031 : ISpergel et al. 



pect there to be significant star formation before this. 

With improving observations, there are many ways to be able to observe these first few 
generations of stars and galaxies. We can now observe high-redshift (z > 6) star formation 
via galactic surveys, which provide statistical properties of high-redshift sources. These 
surveys provide a wealth of information on early star formation; however, they only probe 
the "tip of the iceberg," e.g., a small fraction of the whole population, mainly those galaxies 
that are bright enough to be identified in surveys. 0. 

Another way to observe these high-redshift galaxies is to look for the diffuse back- 
ground light originating from them. Ultra-violet photons produced by these early galaxies 
at z > 6 would be redshifted into the near infrared band (A > 1 /im). Therefore, any 
background above and beyond that from the low-redshift galaxies could be attributable 
to these early galaxies. Observing the Near Infrared Background (NIRB) has the bene- 
fit of probing a population of galaxies as a whole, not just the unusually bright objects 



which can be detected individual 



For this purpose, one can observe the mean NIRB 
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In our previous work ( IFernandez et al.ll2010h . we calculated the angular power spectrum 
of the NIRB fluctuations using an iV-body simulation to trace the large-scale structure and 
the formation of galactic halos, coupled with radiative transfer simulations of reionization. 
Combining this numerical data with the latest analytic modeling of the internal, unresolved 



1 Those exceptionally bright galaxies also might be missed in surveys, becaus e they are rare enough that 
they will not be found within the survey's limited area (jTrenti fc Stiaveliill2008l ) 
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properties of radiation processes such as stellar emission, Lym an-a line, two-photon emission 
free-free and free-bound emission from the formed galaxies ([Fernandez fc Komatsul 12006 
we analyzed the effects of changing the star formation efficiency /*, the escape fraction /, 
and the mass and metallicity of the stars on the NIRB angular power spectrum. 



1 

d— I; 



Our previous simulation was able to resolve dark matter halos down to the minimum 
mass of M min = 2.2 x 10 9 M Q . In the current paper, we extend this analysis by improving 
the mass resolution by more than an order of magnitude to M min = 10 s M Q . We also 
include a new physical effect: the suppression of star formation in small- mass halos due to 
photo-ionization heating (the so-called " Jean s-mass filtering" , which prevents bar yons from 



collapsing into small dark matter halos (e.g., IShapiro et al.lll994t llliev et al.l 120071 )). 



The structure of this paper is as follows. The simulations used for this analysis are 
described in § [21 The fiducial stellar model used is discussed in § [31 and our results are 
presented in § [U This model is then expanded upon in § [51 We compare our results to 
observations in § [6], and conclude in § [71 



2. THE SIMULATION 



In this paper, we use a subset of the simulations presented in llliev et al.l ( 1201 ll ). dis- 
cussed below. These new simulations have a significan t ly hig her mass and spatial resolution 
than the ones we previously used in iFernandez et al.l (120 101 ). The minimum halo mass re- 
solved (with 20 particles or more) is 10 8 M© (vs. ~2 x 10 9 M & before). This much higher 
resolution allows us to explore the effect of the minimum halo mass of the halos that produce 
the ionizing photons on the angular power spectrum. Importantly, we can also study the 
effects of the Jeans-mass filtering of low-mass halos due to the elevated temperature of the 
ionized intergalactic medium (IGM) compared to the neutral phase. This higher temper- 
ature suppresses the future formation of dwarf galaxies and stops the fresh gas infall onto 
already-formed ones. Hence, these new simulations do not simply provide an incremental 
improvement on the resolution, but in fact include a whole new population of so urces which 
are dynamically regulated by the reionization process (see also llliev et al.l 120071 for further 
motivation and details on our model). 

In this work we use two box sizes, 114 h~ l Mpc = 163 Mpc, which we use as our 
fiducial simulations, and 37 h~ l Mpc = 53 Mpc. The latter are computationally much 



2 Hcrc, by the "escape fraction," we refer to a fraction of hydrogen-ionizing photons escaping from the 
galaxies into the intergalactic medium. We ignore dust extinction. 
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Simulation Name 


Box Size 
(Mpc 3 ) 


Minimum 

Halo Mass (M ) 


Suppression 


f-y,large 


f-y, small 




T 


Partial_Supp_HighEff 


163 


10 8 


Yes 


10 


150 


8.3 


0.080 


Partial_Supp_LowEff 


163 


10 8 


Yes 


2 


10 


6.7 


0.058 


No_Supp 


53 


10 8 


No 


0.4 


6 


8.6 


0.078 


Complete_Supp 


53 


10 9 


Yes - complete 


12 





8.3 


0.071 



Table 1: Radiative transfer simulations used in this work. 



cheaper and are used to test a wider range of source properties. The simulation parameters 
are summarized in Table [D Both sets of simulations have the same minimum halo mass 
resolution, M halomin = 1O 8 M . We assume a constant mass-to-light ratio for the halos, 
defined by a parameter / 7 , which gives the number of ionizing photons produced per stellar 
atom which manage to escape from the galaxy into the IGM (jlliev et al.l 120061 ). The low- 
mass (M < 1O 9 M ), suppressible sources and the high-mass (M > 1O 9 M ), unsuppressed 
ones typically have different efficiencies f ltSm aii and f~ h iarge, respectively, reflecting the fact 
that the low-mass halos are expected to have larger populations of Pop III, metal-free stars, 
which are much more efficient emitters of ionizing photons than are Pop II stars. Our 
two fiducial simulations consider the cases of high efficiencies (f^^arge = 10, f^^maii = 150, 
"Partial_Supp_HigfiEff") and low ones (f~ u iarge = 2, f^^maii = 10, "Partial_Supp_LowEff"). 

In addition to our two fiducial cases, we also consider two more complementary cases to 
investigate the more extreme assumptions of either no suppression of low-mass sources by 
reionization (case "No_Supp") or a complete suppression, i.e., assuming that halos with mass 
below 1O 9 M never formed any stars (case "Complete_Supp"). Cases with high efficiencies, 
Partial_Supp_HighEff, No_Supp, and Complete_Supp, reach overlap of ionized regions (by 
design) at approximately the same redshifts, z ov = 8.3, 8.6 and 8.3, respectively, where 
overlap is defined as the time when the mean ionized fraction by mass reaches 99%. The 
corresponding integrated electron scattering optical dept hs are r PS = 0.080, 0.0 78 and 0.071, 



respectively, all in agreement with the WMAP5 results (IKomatsu et al.l l2009h . In contrast, 
the low-efficiency case, Partial_Supp_LowEff, corresponds to a more extended reionization 
history and reaches overlap only by z = 6.7, with correspondingly lower integrated optical 
depth of r = 0.058, slightly below the 2 — a limit from the WMAP 5-year results. 

Throughout this paper, we use the cosmological parameters (f2 m , Qj^, Q h , h) = (0.27, 
0.73, 0.044, 0.7) whi ch are based on the WMAP 5-year results combined with the other 
available constraints (IKomatsu et al.ll2009i ). 
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3. Modeling the Stellar Population 



We model the ionizing sources in the radiative transfer simulations based on the col- 
lapsed halos identified in the structure formation simulations, assuming a constant mass-to- 
light ratio. This is defined by a variable / 7 , which quantifies how many ionizing photons 
escape the parent halo into the IGM over a timestep: 



Here, iVj is the number of ionizing photons produced per stellar atom, / esc is the escape 
fraction, and /* is the star formation efficiency, or the fraction of baryons within the halo 
that is converted into stars. 

For the purposes of modeling the Near Infrared Background, we also need to specify 
the spectral information as well as the spectral energy distribution (SED) of the radiation, 
i.e., we need separate values for JVj and /* above (still satisfying the overall / 7 constraint). 
We shall take Population II stars (with metallicity Z = 1/50 Z Q ) wit h a low escape fraction 
(/esc = 0.1) obeying the following Salpeter initial mass spectrum (jSalpeterl Il955l ) as our 
fiducial case: 

f(m) oc m" 2 - 35 , (2) 

with mass limits of mi = 3 M & and m2 = 150 M , for our entire range of redshifts. The 
star formation efficiency (/*) is then adjusted such that / 7 is held fixed (i.e., the reionization 
constraint is satisfied). The values of /* that are chosen for each value of _f 7 are listed in 
Table [2J 

It is probably unphysical that only these stars will be present at all redshifts, since 
the very first stars were metal-free, with quite different IMF and properties. Nevertheless, 
we make this our fiducial model because the integrated NIRB angular power spectrum sig- 
nal will likely be dominated by the later stages of reionizatiorjfl by which point the metal 
enrichment should be fairly complete. This also represents the case with the highest ampli- 
tude for the angular p ower spectrum among several tested previously (see Figures 7 & 8 of 



Fernandez et al. 



20101 ). We explore other stellar populations in §[5j 



We can then compute the so urce luminosities in the infrare d ban ds following the method - 
ology we presented previously in iFernandez fc Komatsul (120061 ) and Fernandez et al.l (120 lOl ) . 
The total emission has two components, the direct stellar emission; and any light re-processed 



3 Here, we focus on galaxies present during the era of reionization. To observe these galaxies, we assume 
that the low redshift component {z < 6) can be completely removed. In reality, this is a complicated 
observation to make. 
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fl 


Stell 


ar Population 


/* 


fesc 


0.4 


Pop 


II Salpeter 


1.5 


X 


10" 


-3 


0.1 


0.4 


Pop 


III Larson 


1.6 


X 


io- 


-5 


1 


2 


Pop 


II Salpeter 


7.7 


X 


io- 


-3 


0.1 


2 


Pop 


III Larson 


8.0 


X 


io- 


-5 


1 


6 


Pop 


II Salpeter 


2.3 


X 


io- 


-2 


0.1 


6 


Pop 


III Larson 


2.4 


X 


io- 


-4 


1 


10 


Pop 


II Salpeter 


3.8 


X 


io- 


-2 


0.1 


10 


Pop 


III Larson 


4.0 


X 


io- 


-4 


1 


12 


Pop 


II Salpeter 


4.6 


X 


io- 


-2 


0.1 


12 


Pop 


III Larson 


4.8 


X 


io- 


-4 


1 


150 


Pop 


II Salpeter 


5.8 


X 


io- 


-1 


0.1 


150 


Pop 


III Larson 


6.0 


X 


10" 


-3 


1 



Table 2: The properties of the stellar populations. Each population is defined by an escape 
fraction, / csc , which is set to either high (1) or low (0.1). When paired with a stellar mass 
function (Salpeter or Larson) and stellar metallicity (Pop II or Pop III), a value of the star 
formation efficiency, /* is solved for in order to be consistent with reionization. 



by the nebula, including the Lyman— a line, free-free and free-bound emission, and two- 
photon emission. Any ionizing photons that escape the halos are re-processed into nebular 
emission within the IGM. Therefore, as / esc increases, the amount of ionizing photons that 
remain within the halo to produce nebular emission falls, causing the overall luminosity of 
the halo to fall. (These ionizing photons will then produce nebular emission from the IGM, 
but because the IGM is so diffuse, the angular power spe ctrum of the IGM is a nywhere from 
2 to 7 orders of magnitude lower than that of the halo (IFernandez et al.ll2010i ) and thus its 
contribution to the overall NIRB fluctuations can be safely disregarded.) 

These analytical calculations are then combined with the power spectrum of the mass 
density of dark matter halos obtained from our simulations to calculate the 3-dimensional 



power spectrum of the luminosity density, Pl ( k 



, z ) . This 3-D power spectrum is then 



r(z) '■ 

integrated over redshift t o obtain the angular p ower spectrum, which is our key observable 
quantity (see Eq. (37) of IFernandez et all 120101 ): 



C, 



(4tt) 



dz 



H(z)r 2 (z)(l + z) 



-Pl k 



I 



r(z) 



(3) 



where r(z) = c f* dz/H(z) is the comoving distance. The integration range is taken from 
z mm = 6 to z max = 30. The luminosity power spectrum is proportional to the mass density 
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power spectrum of dark matter halos, Pm(/c), as Pi{k) = (L/M) 2 Pm(^)- The constant of 



prop ortionality is the luminosity-to-halo mass ratio given by (see Eq. (9) of Fernandez et al. 
2010h : 

/*^ {**(*) + (1 " /esc) [P f (z) + V*(Z) + P*(z) + V*"(Z)] } • (4) 



M 



a 



Here, l l is the band-averaged l uminosity of a given radiation process % per stellar mass 
calculated from the formulae of Fernandez fc Komatsu ( 20061 ). where i = * (stellar), // 
(free-free), fb (free-bound), 27 (two-photon) and Lya (Lyman-a line). It is important to 
recall two facts: 



1. The power spectrum amplitude is proportional to (L/M) 2 oc f 2 . 

2. The higher the escape fraction / esc is, the smaller the power spectrum becomes. 



Therefore, for a given reionization history (i.e., a given value of / 7 = f* f Pfi r. Nj), the power 
spect rum amplitude is maximized when /* is maximized while / csc is minimized ( Fernandez et al 
2010h . 



4. RESULTS 



4.1. Shape of the Angular Power Spectrum 



In Figure HJ we show the angular power spectra for the following four cases: 1.) A 
minimum dark-matter-halo mass of 10 8 M & with no suppression of the star formation due 
to photo- ionization heating, a minimum halo mass of 10 8 M with 2.) high and 3.) with low 
efficiency with suppression turned on, and 4.) a minimum halo mass of 10 9 M & (in other 
words, complete suppression of low mass sources). Results are presented for the J-band, 
which is assumed to be a rectangular bandpass from 1.1 to 1.4 microns. 



The previous calculations of the angular power spectrum ( IKashlinskv et all2004 : Coorav et al 

20041) yielded a turn-over of l 2 Ci at / ~ 10 3 . Our calculation presented in Fernandez et al. 



( 120101 ) did not confirm this. As we argued in that paper, the previous calculations did not 
take into account the non-linear halo bias, but used the li near bias, underestimat ing the 
power spectrum at high /. See the right panels of Figure 4 of iFernandez et al.l (120101 ) for the 
/c-dependence of the halo bias, b(k), or Figure 6 of Iliev et alT ( 2011 ). 



The effect of non-linear bias will decrease as the minimum mass decreases, and this 
effect can be directly observable. One outstanding question which was left unanswered in 
our previous work is whether a lower minimum halo mass, M min = 10 8 M , can yield a turn 




No Suppression, Min mass=10 M sm 

No Suppression, Min mass=10 9 M sun 

------- Small Halos Suppressed, High Efficiency 

Small Halos Suppressed, Low Efficiency 



Fig. 1. — Angular power spectra of the NIRB fluctuation for Pop-II galaxies (with metallicity 
Z = 1/50 Zq) with a low escape fraction (/ esc = 0.1) and a Salpeter initial mass spectrum. 
The upper solid and dotted lines are for the minimum halo masses of M min = 10 8 M Q 
and 10 9 Mq, respectively, with no suppression (the box size is 53 Mpc). The lower solid 
and dotted lines show the corresponding shot-noise contributions. The upper dashed and 
dot-dashed lines are for the high- and low-efficiency cases (see Table 1), respectively, with 
suppression of star formation in halos with M < 10 9 M & in ionized regions (the box size 
is 163 Mpc). The lower dashed and dot-dashed lines show the corresponding shot-noise 
contributions. The result is shown for the J-band (a rectangular bandpass from 1.1 to 
1.4 /ini). 
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No Suppression, Min mass=10 M sun 

No Suppression, Min mass=10 9 M sun 

------- Small Halos Suppressed, High Efficiency 

Small Halos Suppressed, Low Efficiency 



Fig. 2. — Same as Fig. 1, but for the clustering component of the angular power spectrum 
(Ci minus the shot-noise contribution), normalized to have the same amplitude at I ~ 5000. 
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over in l 2 Ci. We address this question using the high-resolution simulation result presented 
in this paper. In order to see more clearly the effect of the nonlinear halo bias on the 
shape of the angular power spectrum, we show the clustered component (the angular power 
spectrum minus the shot-noise contribution) in Figure |5J To directly compare the shapes of 
these curves, the clustered component of the angular power spectrum for each population is 
normalized to be equal at / pa 5000. 

As expected, we find that the angular power spectrum begins to flatten at high / (around 
i ~ 10 3 ) as the minimum halo mass decreases. This occurs simply because, as the minimum 
halo mass decreases, the non-linear halo bias decreases as well, reducing the small-scale halo 
clustering power, and thus causing the angular power spectrum to flatten. However, at no 
mass above M min = 10 8 M & do we see a turnover of the angular power spectrum. 

The power spectrum is steepest for the (more physically realistic) cases with a par- 
tial suppression of small halos, and for the (extreme) case with a complete suppression of 
M < 10 9 M Q halos. Varying the source efficiencies ( "Partial_Supp_HighEff" versus "Par- 
tial _Supp_LowEff" ) does not change the shape of the power spectra noticeably (but does 
change its amplitude). Therefore, the shape of the angular power spectrum gives information 
on the mass of the halos producing the NIRB fluctuations, with the steepest angular power 
spectrum resulting from populations either with a high minimum halo mass, or where smaller 
galaxies are partially suppressed. Conversely, if low-mass, less biased sources are present and 
not suppressed by reionization, the NIRB angular power spectrum would be much flatter. 

It is important to keep in mind that here, we assume a constant mass-to-light ratio for 
all halos, regardless of mass. In reality, this is a simplified assumption, and it is probable 
that smaller halos will have a smaller mass-to-light ratio. One way to think of this is to 
relate the mass-to-light ratio of the halo to the value of / 7 . If small halos have a smaller 
mass-to-light ratio than larger halos, f 7: i arg e would fall as / 7)Smo // rose, flattening the angular 
power spectrum. We can see the effect of this type of trend if we compare our two cases 
without suppression (the angular power spectrum flattens if we move from f lt i ar ge = 12 and 
fi, small = to f 7i iarge — 0.4 and f y>sm aii — 6)- This could be yet another variable that could 
affect the shape of the angular power spectrum. 



4.2. Amplitude of the Angular Power Spectrum 

Now, we can turn our attention to the amplitude of the angular power spectrum. The 
highest amplitude is produced if only the large sources are active (case "Complete_Supp"), 
i.e., the minimum source mass is 10 9 M & and all smaller halos are completely suppressed. 
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Why does this case with a fewer, massive sources yield a higher NIRB angular power spec- 
trum than does a case when all halos are emitting (case "No_Supp")? The answer is simple: 
when only large halos are active, they are entirely responsible for reionization on their own. 
Consequently, in order to satisfy the reionization constraint, a relatively high ionizing effi- 
ciency of fyjarge = 12 is required. On the other hand, when all halos are emitting (i.e., when 
the minimum source mass is 10 8 M ), the massive source efficiency for the same overlap 
epoch is just / 7l j orse = 0.4. For a fixed initial mass spectrum (hence the number of ionizing 
photons per stellar atom, iVj) and / esc , this translates into 30 times smaller /*, and thus 900 
times smaller contribution to C\. As a result, C\ for this case is dominated by small- mass 
halos. 

Next, we consider the cases with partial suppression of low-mass (< 10 9 M & ) halos due 
to the Jeans-mass filtering. In this instance, small halos will be suppressed in the vicinity 
of large halos; thus, the story is similar to the case with a complete suppression of low-mass 
halos presented above. As some ionizing photons come from low-mass halos, the efficiency 
for large-mass halos is slightly smaller than the case above (f Jy i ar ge = 10 instead of 12; see 
Table 1). This yields the NIRB power spectrum amplitude that is similar to, but slightly 
lower than, the case with M min = 10 9 M & ("Complete_Supp"), and significantly higher than 
the case with M min = 10 s M & without suppression ("No_Supp"). Finally, a test-case scenario, 
in which source efficiencies are so low that the overlap redshift, z ov = 6.7, and the optical 
depth, r = 0.058, are both low, produces a significantly lower angular power spectrum than 
any of the other cases. 



OTHER STELLAR POPULATIONS 



There are still significant uncertainties in the properties of the early stars and galaxies 
- their mass, metallicity, escape fraction, and star formation efficiency - and these can all be 
quite dif ferent from our fiducial case. The dependence on these properties were discussed in 
depth in Fernandez et al.l (120 lOf ). and will only be discussed briefly here. 



There are multiple possibilities available for the stellar population. Along with our 
fiducial case, we will also discuss a case with Pop III (zero-metallicity, Z = 0) stars with a 
high escape fraction (f esc = 1) and a top heavy, Larson initial mass spectrum ( ]Larsonlll998l ): 



f(m) oc m 1 I 1 + 



ni 



-1.35 



(5) 



with m i = 3M Q , m 2 = 500M Q , and m c = 250M o . Our earlier results in Fernandez et al. 
( 120 lOl ) suggested that, while our fiducial case above (Pop II stars with a Salpeter initial mass 
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spectrum and a low escape fraction) represented the limit with the highest amplitude of the 
NIRB angular power spectrum, this new case represents the opposite limit of a low NIRB 
amplitude. Different combinations of escape fraction and star formation efficiency, paired 
with these mass functions and metallicities, will most likely yield amplitudes between these 
two limits. 

We can further vary which populations are found in which type of source - high-mass 
and low-mass halos do not necessarily contain the same stellar populations. The results are 
shown in Figure [31 In all cases, the steepest angular power spectrum results from those cases 
with the largest minimum halo mass of 10 9 M as well as those with a partial suppression, 
and the shallowest from the cases with the lowest minimum halo mass of 10 8 M Q , although 
the detailed shapes and amplitudes vary significantly between different cases. 

The highest NIRB angular power spectrum amplitudes results when the larger sources 
contain the bright Population II stars with low escape fraction. For most cases, the bright- 
est NIRB emission is produced when only the high-mass sources are emitting. In addition, 
the large-mass halos only case "Complete_Supp" is roughly equivalent to the case "Par- 
tial_Supp_HighEff," in which the small halos are suppressed and fyjarge is high. As before, 
this can be explained by the value of f^jarge' angular power spectra with a high amplitude 
have a high value of f lt i a rge- Small halos do not contribute if they are dynamically suppressed. 

Comparison between the top-left (Pop II Salpeter for all halos) and top-right panels (Pop 
II Salpeter for large-mass halos and Pop III Larson for small-mass halos) of Figure 3 would 
show this clearly. As the Pop-III-Larson case has / esc = 1, and its /* is ~ 10 2 times smaller 
than that of the Pop-II-Salpeter case (see Table 2), the small-mass halos are significantly 
dimmer in the top-right panel. The dotted lines are identical in both panels, as these show 
the case with M min = 10 9 M and thus the same Pop II Salpeter. The suppressed cases, 
shown by dashed and dot-dashed lines in both panels, are also very similar. This is because 
small-mass halos are suppressed anyway, and thus the results do not change much by making 
small-mass halos dimmer. On the other hand, the cases with no suppression (the solid lines) 
are very different: in the top- left panel, the small- mass halos dominate the power spectrum 
and thus the power spectrum is much flatter than the dotted line. However, in the top-right 
panel, the small-mass halos are much dimmer and thus the power spectrum is dominated by 
the large-mass halos again, making the shape of the power spectrum as steep as the dotted 
line. (The amplitude is 900 times smaller than the dotted line, as f* t i ar ge for the solid line is 
30 times smaller than that for the dotted line.) 

The one exception is when the massive sources are dimmer, containing Pop III stars 
with a large escape fraction, while low-mass halos contain Pop II stars with a low escape 
fraction. In this case, adding brighter small halos is more than enough to compensate for 
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having "less efficient" large halos, even when these small halos are suppressed. However, 
because of the hierarchical nature of the cosmological structure formation, it is physically 
unlikely to have metal-enriched low-mass halos while having unenriched high-mass halos, 
and thus this case is the least plausible. 



OUR MODEL AND OBSERVATIONS 



There have been multiple recent attempts to measure the fluctuation power of the NIRB 
due to high-redshift, unresolved galaxies. Measurements of the diffuse background light are 
always challenging, as they depend upon careful subtraction of all foreground sources, which 
is complicated. For example, zodiacal light is a main contaminant, and no complete model 
for it is available. When dealing with fluctuations, however, the main uncertainties are due 
to low redshift (z < 6) galaxies. Typically, a map with a certain limiting magnitude is used 
to mask out low redshift galaxies, and another component is used to account for galaxies that 
are too dim to be seen, assuming a certain evolution scenario. In order to obtain accurate 
results, one must be careful that these foreground components are successfully removed while 
still leaving enough of the map to perform accurate fluctuation analysis. 

We have chosen to compare our theoretical predictions with the results of several groups 
who presented measurem ents at a range of frequencies in the near infrared. The first is 
Thompson et al.l (j2007al ). who used the Near Infrared Camera and Multi-Object Spectrome- 
ter (NICMOS) on the Hubble Space Telescopeto measure fluctuations in the NIRB at 1.1 and 
1.6 microns. They removed zodiacal light by dithering t he camera, and removed foreground 
galaxies down to the ACS and NICMOS detection limits. iMatsumoto et al.l (120101 ) measured 
the fluctuations of the NIRB using the AKARI satellite at 2.4, 3.2, and 4.1 microns. The 
camera was rotated over several epochs to remove zodiacal light fluctuations. To remove 
foreground galaxies, they followed a procedure which included masking all pixels over 2 — a, 
repeating the process ten times, and removing foreground sources using DAOFIND. Finally, 
two groups analyzed data from the Inf rared Array Camera (IRAC) on the Spitzer Spac e 
Telescope at 3.6, 4.5, 5. 8 and 8 microns (|Kashlinskv et al.ll2005l l2007bl ; ICooray et al.l 120071 ). 
Kashlinsky et al.l (120051 ) and iKashlinsky et al.l (l2007bl ) removed zodiacal light by taking ob- 
servations 6 months apart in fields ro tated by 180°. Bo t h groups rem oved pixels > 4cr above 
the mean. To remove fainter sources, IKashlinsky et al.l (120051 l2007bt) removed sources fo und 
by SExtractor convolved with the point spread function of IRAC. ICooray et al.l ( 120071 ). on 
the other hand, masked the image to various magnitude limits and used known galaxies in 
the ACS catalog. 



In Figure HI we compare our fiducial calculations (with Pop II stars; a Salpeter initial 
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No Suppression, Min mass=10 M sun 

No Suppression, Min mass=10 9 M sun 

------- Small Halos Suppressed, High Efficiency 

Small Halos Suppressed, Low Efficiency 



Fig. 3. — Same as Fig. 1, but for various stellar populations. Top left panel: the same 
population as Fig. 1 (Pop II Salpeter for all halos). Top right panel: Pop II Salpeter for 
large-mass (> 10 9 M ) halos; Pop III Larson for small-mass (< 10 9 M ) halos. Bottom left 
panel: Pop III Larson for large-mass halos; Pop II Salpeter for small-mass halos. Bottom 
right panel: Pop III Larson for all halos. 
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Fig. 4. — Our fiducial model (Pop II stars with a Salpeter initial mass spectrum and f e . 



bands, compared 


to the 


(Thompson et al. 


2007a. 



top-left panel), 2.4 m icrons from AKARI (Matsumoto et al. 2010l top- r ight panel), and 3.6 
microns from IRAC kashlinsky et al.ll2005l l2007c! : ICoorav et al.ll2007l bottom-left panel). 
For our models, the top line in each case is without the shot-noise subtracted, and the 
bottom line is with our p rediction for the sho t-noise subtracted. Note that in some cases 
(ie, the measurements for ICoorav et al.l ( 120071 )). a higher model of shot-noise than ours is 
subtracted from the observations (namely, where the shot-noise is fit to the angular power 
spectrum at high multipoles, and the data may contain contributions from both high-;? and 
low-z galaxies). 
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mass spectrum; and f PSC = 0.1) to the measur ed angular power spectra at 1.6 microns from 



the N ICMOS camera ([Thompson et al.ll2007af). 2.4 microns from AKARI (IMatsumoto et al. 



20101), and 3.6 microns from IRAC ( IKashlinskv et al.ll2005l . l2007ri ICoorav et al.ll2007f ). In all 



cases, the measurements are above our calculations. These measurements should probably 
be taken as upper limits on the contributions from galaxies in z > 6, as there may be 
some lower redshift components or other foreground contamination in these measurements. 
We also note that our fiducial model presented here has the largest predicted amplitude 
of Cf, thus, if we were to change the population within our range of models satisfying 
the reionization constraint, the predicted amplitude of the angular power spectrum would 
decrease. In other words, these upper limits do not yet rule out any of our models. We also 
note that current observational limits are many orders of magnitude above that we would 
expect from our lowest magnitude of the angular power spectrum. Many of our models are 
far beneath the sensitivity of curr ent observations. T he upcoming experiment CIBER will 
push the sensitivity limit further ( ICoorav et al.l 120091 ). and may begin to be able to discern 
between various high-redshift populations. Still, it is possible that many o f our models will 



need e xtraordinarily sensitive instruments to be detected (see figure 18 of Fernandez et al. 



7. CONCLUSIONS 

The NIRB mean and fluctuation signals can be a very powerful prob e of the high-redshift 



star formation (z > 6). Build ing upon our previous work on this subject ( Fernandez &: Komatsu 



20061 ; Fernandez et al.ll2010l ). we have extended the calculation of the angular power spec- 
trum of the NIRB fluctuation by improving the mass resolution of A^-body simulations by an 
order of magnitude. With the simulation now resolving the dark matter halo mass down to 
M = 10 8 Mq, we have confirmed our previous finding: due to the non-linear halo biasing, the 
shape of l 2 Ci of the NIRB fluctuation does not exhibit a turn over. However, we do observe 
a flattening of the shape as we lower the minimum halo mass (because the non-linear bias 
decreases as masses go down), and thus the shape of l 2 C[ can be used to infer typical masses 
of dark matter halos hosting galaxies in a high redshift universe. 

We have gone beyond simply increasing the mass resolution of the simulation. For 
the first time, we consider cases with a radiative feedback suppressing the star formation 
in low-mass galaxies with M < 10 9 M , due to the Jeans-mass filtering in the ionized 
and heated IGM. We find that, when low- mass sources are partially suppressed by photo- 
ionization heating, the predicted angular power spectrum becomes similar to the one with 
a complete suppression of halos below M min = 10 9 M , yielding a steep power spectrum at 
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high I. Therefore, the shape of the angular power spectrum at high I can directly provide 
information on the typical mass of sources responsible for completing reionization. 

The amplitude of the angular power spectrum is less robust, as it depends on a number 
of parameters such as the stellar initial mass spectrum, metallicity of stars, star formation 
efficiency, and escape fraction. However, one robust feature is that it is largely determined by 
the properties of the halo populations at late stages of reionization. Therefore, the angular 
power spectrum is higher if the stars produce more ionizing photons that do not escape from 
the halo. (Therefore, maximizing /* and iVj while minimizing / osc within the limitations of 
reionization.) 

Finally, we find that our predictions, all of which are tuned to satisfy the reionization 
constraints, are below the current measurements. Given that these measurements should 
probably be taken as upper limits on the contributions from galaxies in z > 6, we conclude 
that our calculations are consistent with the current measurements. 
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